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We present results for the charmonium-nucleon potential V c e-n{t) from quenched lattice QCD, 
which is calculated from the equal-time Bethe-Salpeter amplitude through the effective Schrodinger 
equation. Detailed information of the low-energy interaction between the nucleon and charmonia 
(r/ c and J/tp) is indispensable for exploring the formation of charmonium bound to nuclei. Our 
simulations are performed at a lattice cutoff of 1/a ~ 2.1 GeV in a spatial volume of (3 fm) 3 with 
the nonperturbatively 0(a) improved Wilson action for the light quarks and a relativistic heavy 
quark action for the charm quark. We have found that the potential V C c-n(t) for either the r\ c and 
J/tp states is weakly attractive at short distances and exponentially screened at large distances. The 
spin averaged J/ip-N potential is slightly more attractive than that of the rj c -N case. 



Heavy quarkonium states such as charmonium (cc) 
states do not share the same quark flavor with the nu- 
cleon (N). This suggests that the heavy quarkonium- 
nucleon interaction is mainly induced by the genuine 
QCD effect of multi-gluon exchange [THS]- Therefore the 
cc-N system is ideal to study the effect of multi-gluon 
exchange between hadrons. As an analog of the van der 
Waals force, the simple two-gluon exchange contribution 
gives a weakly attractive, but long-ranged interaction be- 
tween the heavy quarkonium state and the nucleon j^l [5] . 
However, the validity of calculations based on a pertur- 
bative theory is questionable for QCD where the strong 
interaction influences the long distance region. 

The cc-N scattering at low energy has been studied 
from first principles of QCD. The s-wave J/ip-N scatter- 
ing length is about 0.1 fm by using QCD sum rules [6] 
and 0.71 ± 0.48 fm (0.70 ± 0.66 fm for r/ c -N) by lattice 
QCD [7], while it is estimated as large as 0.25 fm from the 
gluonic van der Waals interaction [2] . All studies suggest 
that the cc-N interaction is weakly attractive. This indi- 
cates that the formation of charmonium bound to nuclei 
is enhanced. In 1991, Brodsky et al. had argued that 
the cc-nucleus (A) bound system may be realized for the 
mass number A > 3 if the attraction between the charmo- 
nium and the nucleon is sufficiently strong ^ [5] . There- 
fore, precise information on the cc-N potential V C c-n(t) 
is indispensable for exploring nuclear-bound charmonium 
states like the 7y c - 3 He or J/-0- 3 He bound state in few body 
calculations [5]. 

We recall the recent great success of the N-N potential 
from lattice QCD [TU]. In this new approach, the poten- 
tial between hadrons can be calculated from the equal- 
time Bethe-Salpeter (BS) amplitude through an effective 
Schrodinger equation. Thus, direct measurement of the 
cc-N potential is now feasible by using lattice QCD. It 
should be important to give a firm theoretical prediction 
about the nuclear-bound charmonium, which is possibly 
investigated by experiments at J-PARC and FAIR/GSI. 

The method utilized here to calculate the hadron- 
hadron potential in lattice QCD is based on the same 



idea originally applied for the N-N potential [101 fTT] , 
We first calculate the equal-time BS amplitude of two 
local operators (hadrons hi and h 2 ) separated by given 
spatial distances (r = |r|) from the following four-point 
correlation function: 

G h ^(r,t;t 2 ,t 1 ) 

x x',y' 

where r is the relative coordinate of two hadrons at sink 
position (i). Each hadron operator at source positions 
(ti and t 2 ) is separately projected onto a zero-momentum 
state by a summation over all spatial coordinates x' and 
y'. To avoid the Fierz rearrangement of two-hadron op- 
erators, it is better to set t 2 ^ t%. Without loss of gen- 
erality, we choose t 2 = ti + 1 = t STC hereafter. Suppose 
that \t — t STC \ 3> 1 is satisfied, the correlation function 
asymptotically behaves as 

G hl ~ h2 (r, i; i src ) oc 4> hl -h 2 (r) e - E ^ (1) 

where the r-dependent amplitude, which is defined as 

K-M=^{0\hi(x)h 2 (x + T)\h 1 h2',E hl . Ka ) (2) 

x 

with the total energy E^.^ for the ground state of the 
two-particle hi-h 2 state, corresponds to a part of the 
equal-time BS amplitude and is called the BS wave func- 
tion [T2 EI]- After an appropriate projection with re- 
spect to discrete rotation 

<^» = ^ E ^-hdK-'r), (3) 

where 71 represents 24 elements of the cubic group Oh, 
one can get the BS wave function projected in the Af 
representation, which corresponds to the s-wave in con- 
tinuum theory at low energy. Once the BS wave func- 
A + 

tions 4> h % _ h (r) are calculated in lattice simulations, the 
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TABLE I: Simulation parameters in this study. The Sommer 
parameter ro = 0.5 fm is used to fix the scale [141 115] . 
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Statistics 


6.0 


0.0931 


2.12 


16 a x 
32 3 x 


48 
48 


1.5 
3.0 


200 
602 



hadron-hadron "effective" central potential with the en- 
ergy eigenvalue E of the stationary Schrodinger equation, 
can be obtained by 

1 V 2 <^\ (r) 

where fj, is the reduced mass of the hi~h 2 system and V 2 is 
defined by the discrete Laplacian with nearest-neighbor 
points. Although the energy eigenvalue E is supposed 
to be the energy difference between the total energy of 
two hadrons {Eh y -h 3 ) an d the sum of the rest mass of 
individual hadrons (M^ + M/j 2 ), we instead determine 
E with the condition oi ]im r - yoo {Vh 1 -h 2 ( r ) — E} = [15] . 
More details of this method can be found in Ref. [11] . 

Let us first consider the low energy rj c -N interac- 
tion, which does not possess a spin-dependent part. We 
use the conventional interpolating operators, h\(x) — 
£abc(u a (x)C"f5db(x))u c (x) for the nucleon and Ii2(y) — 
Ca{y)l5C a {y) for the r\ c state, where a, b and c are color 
indices, and C = 7472 is the charge conjugation matrix. 
We have performed quenched lattice QCD simulations on 
two different lattice sizes, L 3 x T = 32 3 x 48 and 16 3 x 48, 
with the single plaquette gauge action at /? = 6/g 2 = 6.0, 
which corresponds to a lattice cutoff of a -1 « 2.1 GeV 
according to the Sommer scale [TH [TS] . Our main results 
are obtained from the data taken on the larger lattice 
(La « 3.0 fm). A supplementary data with a smaller 
lattice size (La ~ 1.5 fm) are used for a test of finite-size 
effects. The number of statistics is O(600) for L = 32 and 
O(200) for L = 16, respectively. The simulation parame- 
ters and the number of sampled gauge configurations are 
summarized in Tabled] 

We use non-perturbatively 0(a) improved Wilson 
fermions for light quarks (q) 16] and a relativistic heavy 
quark (RHQ) action for the charm quark (Q) [17] , 
The RHQ action is a variant of the Fermilab approach 
[18] . which can remove large discretization errors for 
heavy quarks. The hopping parameter is chosen to be 
K q = 0.1342, 0.1339, 0.1333, which correspond to = 
0.64, 0.72, 0.87 GeV (M N = 1.43, 1.52, 1.70 GeV), and 
kq = 0.1019 which is reserved for the charm-quark mass 
(M Vc = 2.92 GeV and M J/4 , = 3.00 GeV) 19 . Each 
hadron mass is obtained by fitting the corresponding 
two-point correlation function with a single exponential 
form. We calculate quark propagators with wall sources, 
which are located at t STC — 5 for the light quarks and 



at t src = 4 for the charm quark, with Coulomb gauge 
fixing. It is worth mentioning that Dirichlet boundary 
conditions are imposed for quarks in the time direction 
in order to avoid wrap-round effects, which are very cum- 
bersome in systems of more than two hadrons [20] , In 
addition, the ground state dominance in four-point func- 
tions is checked by an effective mass plot of total energies 
of the cc-N system. 

The left panel of Fig. [I] shows a typical result of the 
projected BS wave function at the smallest quark mass, 
which is evaluated by a weighted average of data in 
the time-slice range of 16 < t — t src < 35. The wave 
functions are normalized to unity at a reference point 
r = (16, 16, 16), which is supposed to be outside of the 
interaction region. As shown in Fig. [T] the wave func- 
tion is enhanced from unity near the origin so that the 
low-energy rj c -N interaction is certainly attractive. This 
attractive interaction, however, is not strong enough to 
form a bound state as is evident from this figure, where 
the wave function is not localized, but extends to large 
distances. 

In the right panel of Fig. [T] we show the effective cen- 
tral rj c -N potential, which is evaluated by the wave func- 
tion through Eq. Q with measured values of E and fi. 
As is expected, the r] c -N potential clearly exhibits an en- 
tirely attractive interaction between the charmonium and 
the nucleon without any repulsion at either short or large 
distances. The short range attraction is deemed to be a 
result of the absence of Pauli blocking, that is a relevant 
feature in this particular system of the heavy quarkonium 
and the light hadron. The interaction is exponentially 
screened in the long distance region r > 1 fm [25] . This is 
consistent with the expected behavior of the color van der 
Waals force in QCD, where the strong confining nature of 
color electric fields must emerge [5, 21 . The exponential- 
type damping in the color van der Waals force is hardly 
introduced by any perturbative arguments. 

In detail, a long-range screening of the color van der 
Waals force is confirmed by the following analysis. We 
have tried to fit data with two types of fitting functions: 
i) exponential type functions — exp(— r m )/r™, which in- 
clude the Yukawa form (m = 1 and n = 1), and ii) in- 
verse power law functions —l/r™, where n and m are not 
restricted to be integers. The former case can easily ac- 
commodate a good fit with a small x 2 /ndf value, while 
in the latter case we cannot get any reasonable fit. For 
example, the functional forms — exp(— r)/r and —l/r 7 
give x 2 / n df — 2.5 and 34.3 , respectively. It is clear 
that the long range force induced by a normal "van der 
Waals" type interaction based on two-gluon exchange [5] 
is non-perturbatively screened. 

If we adopt the Yukawa form —-fe~ ar /r to fit our data 
of V c£ -N(r), we obtain 7 ~ 0.1 and a ~ 0.6 GeV. These 
values should be compared with the phenomenological cc- 
N potential adopted in Refs. .1], where the parameters 
(7 = 0.6, a = 0.6 GeV) are barely fixed by a Pomeron 
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FIG. 1: The wave function (left) and the effective central potential (right) in the s-wave r) c -N system for M v = 0.64 GeV as 
a typical example. In the right panel we fit a Yukawa potential (solid line) and compare with the phenomenological potential 
(dashed) adopted in Ref. [JJ. 



exchange model. The strength of the Yukawa potential 
7 is six times smaller than the phenomenological value, 
while the Yukawa screening parameter a obtained from 
our data is comparable. The cc-N potential obtained 
from lattice QCD is rather weak. 

We next show both finite-size and quark-mass depen- 
dence of the rj c -N potential in Fig. [2j Firstly, as shown 
in the left panel of Fig. [2] there is no significant differ- 
ence between potentials computed from lattices with two 
different spatial sizes (La ps 3.0 and 1.5 fm). This obser- 
vation is simply due to the fact that the rj c -N potential 
is quickly screened to zero and turns out to be some- 
how short ranged. In principle, the short range part of 
the potential, which is represented by ultraviolet physics, 
should be insensitive to the spatial extent associated with 
an infrared cutoff. As a result, it is assured that the larger 
lattice size is large enough to study the r) c -N system. 

No large quark-mass dependence is also observed in 
the right panel of Fig. [2j This is a non-trivial feature 
since there is an explicit dependence on the reduced mass 
fi in the definition of the effective central potential Q. 
However, if one recalls that the cc-N interaction is mainly 
governed by multi-gluon exchange, the resulting potential 
is expected to be less sensitive to the reduced mass of the 
considered system ignoring the internal structures of the 
t] c and nucleon states. 

If one takes a closer look at the inset of Figure [2j it 
is found that the nature of the attractive interaction in 
the r] c -N system tends to get slightly weaker as the light 
quark mass decreases. Does this mean that the strength 
of the rj c -N potential at the physical point becomes much 
weaker than what we measured at the quark mass sim- 
ulated in this study? The answer to this question is not 
simple. It is worth to remember that the ordinary van der 



Waals interaction is sensitive to the size of the charge dis- 
tribution, which is associated to the dipole size. Larger 
dipole size yields stronger interaction. We may expect 
that the size of the nucleon becomes large as the light 
quark mass decreases. However, the very mild but oppo- 
site quark-mass dependence observed here does not ac- 
commodate this expectation properly. 

Recent detailed studies of nucleon form factors tell us 
that the root mean-square (rms) radius of the nucleon, 
which is a typical size of the nucleon, shows rather mild 
quark-mass dependence and its value is much smaller 
than the experimental value up to at ~ 0.3 GeV 
(for example, see [23]). At the chiral limit in baryon 
chiral perturbation theory the rms radius is expected to 
diverge logarithmically 24J. This implies that the size 
of the nucleon increases drastically in the vicinity of the 
physical point. It may be phenomenologically regarded 
as the "pion-cloud" effect. We notice here that our simu- 
lations are performed in the quenched approximation and 
at rather heavy quark masses. We speculate that the cc- 
N potential from dynamical simulations would become 
more strongly attractive in the vicinity of the physical 
point, where the size of the nucleon is much larger than 
at the quark mass simulated in this study. 

We also have calculated the J/ip-N potential. It should 
be noted here that different total spin states are al- 
lowed as spin-1/2 and spin-3/2 states in the J/ip-N sys- 
tem. This fact introduces a little complexity regard- 
ing the spin-dependence on the J/ip-N potential. The 
interpolating operator of the J/ip state is defined as 
h2,i(y) = Ca(y)liC a (y) , which carries the specific spin 
polarization direction. Therefore, the four-point correla- 
tion function for the J/ip-N state becomes a matrix form 
with spatial Lorentz indices, G^ 1 "' 12 = (hih2,i(hih2 ! j)^)- 
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FIG. 2: The volume dependence (left) and the quark-mass dependence (right) of the rj c -N potential. 
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FIG. 3: The central and spin-independent part of the J/ip- 
N potential at M n = 640 MeV. The r/ c -N potential is also 
included for comparison. 



It can be expressed by an orthogonal sum of spin-1/2 and 
spin-3/2 components f7]: G^ h2 = G^ 2 V lj 2 + G 3 / 2 VfJ 2 
where proper spin projection operators for the spin-1/2 
and spin-3/2 contributions are given by — 
and VfJ 2 = Sij — ^Ji^j in the center of mass frame [2"2"] . 
Then, each spin part can be projected out as G 1 / 2 = 

v.X^v h ' and G3/2 = lEw^S 72 ^ where 

the indices i and j are also summed over all spatial di- 
rections. 

As a result, we can obtain the BS wave function and 
the resulting J/ip-N potential for each spin channel. Al- 
though the lower-spin state (spin-1/2) is not free from the 
contamination of the r) c -N state through channel mix- 
ings [7] , we simply consider the spin averaged four-point 



correlation function for the J/ip-N system as 

G J Jt N = \G 1/2 + \g*' 2 = \ £ G J ^- N . (5) 

i 

This procedure may provide only a spin-independent part 
of the J/ip-N potential through the same analysis applied 
to the ^e-A^ system. 

We show our result of the J/tp-N potential in Fig. [3] 
where the rj c - N potential is included for comparison. The 
J/ip-N potential shows short-range attraction. Similar 
to what we discussed in the rj c -N case it does not have 
a normal "van der Waals type" — 1/r" behavior. There 
is no qualitative difference between the ri c -N and J/ip- 
N potentials. Quantitatively, the attractive interaction 
observed in the J/ip-N potential is rather stronger than 
that of the r] c -N system as shown in Fig. [3j though it 
is still not strong enough to form a bound state in the 
J/t/j-N system. 

What is a possible origin of the stronger attraction ap- 
pearing in the J/ip-N system? As was discussed previ- 
ously, the attractive interaction in the cc-N system tends 
to be slightly stronger as the light quark mass increases. 
It should be recalled that the reduced mass of the cc-N 
system is also changed through a variation of the light 
quark mass. Supposed that the strength of the cc-N po- 
tential depends simply on the reduced mass of the cc-N 
system, a mass difference of the rj c and J/ip state may 
account for the difference between the r] c -N and J/tp- 
N potentials. However, this is not the case. As shown 
in the inset of Figure [2j the r\ c - N potential gets deeper, 
when the reduced mass of the r\ c -N system is increased by 
about 10% through a variation of the light quark mass. 
On the other hand, although the reduced mass receives 
only about 1% gains in the J/tp-N system relative to the 
rj c -N system, the difference between the rj c -N and J/ip- 
N potentials shown in Fig. [3] is much bigger than the 
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reduced mass dependence observed in Fig. [2j This may 
indicate that the stronger interaction in the J/ip-N sys- 
tem than in the r] c -N system is caused by some dynamics 
associated with the structure of quarkonia. 

We again recall that the ordinary van der Waals inter- 
action is sensitive to the size of the charge distribution. 
Therefore, what we observed here is intuitively accounted 
for by the simple speculation that the size of the J/tp 
state is larger than that of the r] c state. Nevertheless, 
similar to the r) c -N system in the J/ip-N no appreciable 
finite-size effects are observed. 

We have studied the charmonium-nucleon potential 
V cS -n(t) in quenched lattice QCD. Potentials between 
charmonia (r/ c and J/tp) and the nucleon are calculated 
from the equal-time BS amplitude through the effective 
Schrodinger equation. We have found that the central 
and spin-independent potential V C c-n{t) m both the rj c - 
N and J/ip-N systems is weakly attractive at short dis- 
tances and exponentially screened at large distances. It is 
observed that both potentials have no appreciable finite- 
size dependence and no significantly large quark-mass de- 
pendence within the pion mass region 640MeV < M n < 
870MeV. The J/ip-N system is slightly stronger attrac- 
tive than the rj c -N system. This should be accounted for 
by the dynamics associated with the structure of quarko- 
nia. 

In order to make a reliable prediction about the 
nuclear-bound charmonium, an important step in the fu- 
ture is clearly an extension to dynamical lattice QCD 
simulations in the lighter quark mass region. Such plan- 
ning is now underway. Once we obtain a realistic poten- 
tial, we will proceed exploring the nuclear-bound charmo- 
nium state with theoretical inputs of the charmonium- 
nucleon potential by using the exact few-body calcula- 
tion. 
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